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Abstract 

We investigate one-loop renormalization group evolutions of extended sectors of 
Yukawa type couplings. It is shown that Landau Poles which usually provide nec¬ 
essary low energy upper bounds that saturate quickly with increasing initial value 
conditions, lead in some cases to the opposite behaviour: some of the low energy 
couplings decrease and become vanishingly small for increasingly large initial con¬ 
ditions! We write down the general criteria for this to happen in typical situations, 
highlighting a concept of repulsive quasi-fixed points, and illustrate the case both 
within a two-Yukawa toy model as well as in the minimal supersymmetric stan¬ 
dard model with R-parity violation. In the latter case we consider the theoretical 
upper bounds on the various couplings, identifying regimes where ^'kkk^ ^ski 
are dynamically suppressed due to the Landau Pole. We stress the importance of 
considering a large number of couplings simultaneously. This leads altogether to 
a phenomenologically interesting seesaw effect in the magnitudes of the various R- 
parity violating couplings, complementing and in some cases improving the existing 
limits. 
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1 Introduction 


The Infrared Quasi Fixed Point (IRQFP) structure in the Yukawa sector, |^, has played 
an important role in singling out the supersymmetric extension of the standard model 
in the small tan (3 regime, as a natural framework to accommodate a top quark mass of 
~ 175 GeV [@|, as opposed to the standard model itself or to alternatives such as the 
gauged Nambu-Jona-Lasinio ^ model, which predicted altogether too heavy a top quark 
( > 200 GeV). The experimental (almost) negative search of Higgs particle pushes it’s 
mass above 114 GeV in the standard model scenario @, |^, or above 89 GeV |^, in the 
Minimal Supersymmetric Standard Model (MSSM), ruling out the small tan/? IRQFP 
regime. This motivates naturally the study of supersymmetric scenarios where non zero 
bottom and r quark Yukawa couplings are of the same order of magnitude as that of the 
top 10. By the same token, one can also consider further extensions such as the next 
to minimal supersymmetric standard model ((M+1)SSM) P, or the MSSM without R- 
parity 0 , leading to an increasing number of Yukawa type couplings. 

The aim of the present paper is to study the generic behaviour of the renormalization 
group evolutions of the Yukawa couplings in the presence of a large number of these cou¬ 
plings, and in particular the effect of Landau Poles on the perturbative bounds. The main 
result is that for increasingly large values of the initial conditions, some of the couplings 
can be suppressed (formally to zero!) at low energies, a behaviour which differs drastically 
from the usually expected one 0. This phenomenon is related to the appearance of a rich 
structure of Infrared Quasi Fixed Points, some of which can have an unconventional re¬ 
pulsive character. In practice this means that the close numerical connection between the 
Landau Pole and the perturbativity bound, usually observed for the top quark Yukawa 
coupling for instance, does not prevail for other couplings in more general conhgurations. 


A natural place where this new mechanism occurs is the MSSM with R-parity violation 
(/?p-MSSM). Here one has 45 Yukawa couplings as new free parameters to be added to 
those of the MSSM. In the absence of theoretically compelling reasons to put them to zero 
(apart from the elegant [|^, yet ad hoc, R-parity!) the new parameters have to be allowed 
for and studied phenomenologically for their own sake. They can influence many physical 
processes and have already been experimentally constrained this way . An important 
motivation for constraining R-parity violating scenarios is the bearing it can have on the 
nature of the supersymmetric dark matter candidate in the Universe, and indirectly on 
the origin of supersymmetry breaking. While conserved R-parity stabilizes the lightest 
supersymmetric particle (LSP) and tends to favour the lightest neutralino, rather than 
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the gravitino, as the LSP (a natural configuration in gravity mediated supersymmetry 
breaking scenarios), even an extremely small violation of R-parity becomes a ’’Damocles 
sword” over the neutralino LSP as a serious candidate for cold dark matter. Conversely, 
a moderate amount of R-parity violation would favour the gravitino as a potential dark 
matter candidate [O , and to some extent favour scenarios with low scale supersymmetry 
breaking. 


The mo del-independent mechanism we elucidate in this paper can have interesting impli¬ 
cations on the study of R-parity violation. The rest of the paper will be thus divided into 
two main parts -the first part, corresponding to sections 2 and 3 and appendices A and 
B, is devoted to a detailed discussion of the generic solutions of the one-loop renormaliza¬ 
tion group equations for an arbitrary number of Yukawa couplings, and in particular the 
dynamical attraction to null couplings due to the presence of Landau Poles -the second 
part, i.e. section 4 and appendix C, is an illustration of these features in the context 
of |?p-MSSM. There we study the effect of increasing the number of active couplings by 
considering subsets with 6, 9 and 13 non-zero couplings. We compare our results with 
previous theoretical bounds as well as with experimental limits. In appendix C we write 
explicitly the one-loop renormalization group equations governing the evolution of a large 
set of couplings. In section 5 we conclude and summarize our results. [The reader inter¬ 
ested mainly in phenomenological implications on R-parity violation can go directly to 
section 4 which is fairly self-contained.] 

2 General Formulation 

2.1 Integrated Forms for the Running Yukawa Couplings 

Let us consider first a general gauge theory with an arbitrary number of gauge group 
factors and of fermion multiplets coupled to some scalar fields. Through this paper we 
will always assume no flavour violating fermion-fermion-scalar Yukawa couplings. The 
renormalization group equations (RGE) for the gauge couplings, and [with some further 
simplifying assumptions about the scalar self couplings] those for the Yukawa couplings, 
take the following form at the one-loop level |]T^, 


A 

dt 


a,; 


d ~ . 



i I 
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( 1 ) 

( 2 ) 





where t = Log[MQjjrp /denotes the scale evolution parameter, and we dehne 




( 3 ) 


where hk,gi denote respectively the Yukawa and gauge couplings. The Cki and aui are 
constant coefficients depending on the model. 

The general solution for such a system (valid for an arbitrary number of Yukawa couplings 
labeled by k) reads [pT 


ai{t) 

n(t) 


Qi(0) 

1 + 6idi(0)f 

Yk{^)uk(t) 

1 + afcfelfc(O) jQUk{t')dt' 


( 4 ) 


where the auxiliary functions Uk are given by 


Uk{t) = 


Ek{t) 


}}k Jo 


and the functions Ek{t) = ea:p[/Q J2iCkiOii{t')dt'] read explicitly 


( 5 ) 


^ <^ki 

Ek{t) = 11(1 + 6idi(0)t) '’i (6) 

i=l 

In we called these solutions “integrated forms” to stress the fact that they are not 
explicit but rather iterative. As it turned out, and will be illustrated once more in this 
paper, they allow, even in this form, a generic extraction of several interesting features, 
which could be hardly pinned down from mere numerics, in particular in relation with 
the infrared quasi hxed points. In many physically interesting cases such as the non 
supersymmetric standard model, or the MSSM or the NMSSM, the RGB’s for the Yukawa 
couplings take indeed the form of Eq. (||) and the above solutions are directly applicable. 
This is no more strictly true in the general case of ,^p-MSSM, nonetheless, we will see 
that these solutions still prove very efficient in describing this case too. 


2.2 Asymptotic behaviour and quasi fixed points 

Avoiding Landau poles in the Yukawa system leads to consistency upper bounds on the 
values of the running Yukawa couplings at some low energy scale. These bounds can be 
most straightforwardly determined by looking at the asymptotic behaviour of the running 
couplings when some, or all, of the initial conditions are taken inhnitely large. We will 
come back at length to the meaning of these bounds in the subsequent sections. Here 
we derive for further reference the complete procedure which allows a non ambiguous 
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determination of the asymptotic behaviour, recalling and extending the point made in 
ref. Ol. It is useful to dehne 


n(o) 

^.( 0 ) 


(7) 


and study the asymptotic behaviour by increasing the hi( 0 )’s while keeping the above 
scaling ratios hxed and hnite ( 7 ^ 0). One can thus carry out the discussion in terms of the 
set of hnite r’s and one single initial condition parameter Y ( 0 ) = which is allowed to 
become inhnite. In this limit the values taken by the Yk{t) at a low energy scale t are the 
so-called infrared effective hxed points or Infrared Quasi Fixed Points (IRQFP)|5 From 
Eq.(l) one expects the u’s to behave like 


,00 _ 


QFP 


(Oi) 


for asymptotically large Y^, where we will refer to the p^’s as “asymptotic powers”. Here 
is an integrated form which is Y^ independent but may or may not depend on the 
scaling parameters r. Now it is important to stress that Eq.(|D is quite general and does 
not mean|^ that we neglect the 1 in the denominator of Eq. (|^) due to increasingly large 
F^(0)’s. It is indeed easy to see from Eqs.(|, H) that the IRQFP are controlled by the 
various asymptotic powers in the following way, 




akk IoUk^^{t')dt' 
1 + akk lo 


a Pk < I 

if Pfc = 1 


0 


if pk>l 


(9) 


Only when < 1 is it justihed to drop 1 in the denominator. 


In this case the IRQFP’s 


take the usual form as it was shown for the MSSM |^. The two other cases in Eq.(^ 
occur for larger Yukawa sectors as it was found for the (M-l-l)SSM in |]^ and for the ^p- 
MSSM in M]. In particular, the very unusual conhguration where some IRQFP’s vanish 


is worth attention. As was stressed in |T^, the prior determination of the asymptotic 


^ Let us mention briefly here a different definition adopted in [^, where one requires ^yfc('r)|r=t — 0 
for all k, at a given low scale t and solves the resulting linear system of equations in the positive Yk{tys 
from Eqs.(|^). There is, however, no guarantee in general that this linear system has positive solutions; 
one then has to change accordingly the amount of “zeroing” of ^Ffe(T) or/and the scale t. Our definition 
is more general in that the slowing down of the running is automatically such that the IRQFP’s always 
exist 

^With that respect an unfortunate erroneous statement slipped in refs. @ 0 - 
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powers is crucial. Straightforward inspection shows that they must be solutions of the 
following equation, 

V = M-{l-V)e{l-V)5 (10) 

Here P is a column vector of all the asymptotic powers pk'-, (1 — 'P)9{1 — P) is a shorthand 
for the column vector with components (1 — pk)0{l — Pk)5k where 9 is the Heaviside 
function, = 1 or 0 respectively for ^^(O) inhnite or hnite; M. is dehned by 




( 11 ) 


^33 


Note that due to the 9 function in Eq.(pi)D this system is linear in the p^’s only in patches. 
It can thus allow simultaneously for more than one set of solutions which should corre¬ 
spond to alternated conhgurations given by Eq.(|^). We will study at length the meaning 
of such multiple solutions in the following sections, highlighting the (new) phenomenon 
of repulsive IRQFP. 


2.3 An example: the two-Yukawa system 

It is instructive to illustrate the above as well as the forthcoming features of the paper in 
the simplest case of two Yukawa couplings (the “top/bottom” system): 


^Yi(t) = Hi(t)(^ cugiit) - aiiYi(t) - ai2Y2{t)) 

—Y2{t) = Y2{t)C^C2igi{t) — a2iYi{t) — a22Y2{t)) ( 12 ) 

where generically au = 021 = |a|,aii = 022 = |&| as in the MSSM (or in the Standard 
Model (SM)). We will stick to this conhguration for simplicity although our discussion 
applies more generally. In this case Eq. (0) reads 




(1 -pi) 9[l-pi] 
(1 -P 2 ) 9[1-P2] 


(13) 


The critical issue here is whether ||| < or > 1. As can be seen from Appendix B.l, in the 
hrst case Eq.(|T^) has a unique solution and thus one single Infrared Quasi Fixed Point 
given by the hrst equation in (^. This case corresponds to realistic models like the SM 
and the MSSM (where ||| = 1/3 and 1/6 respectively). In the “toy model” case ||| > 1, 
comparing Eqs.([T3|) and (EH (i.e. a = 6 = 0 and f3 = = |a|/|6|), one hnds (see 

Eqs. (p.l.8|, |B.1.9|) ) that there are three different solutions to Eq. (p!3[) corresponding to 
three possible IRQFP’s with the appropriate conhgurations in Eq.(|^). Which one of these 
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solutions is dynamically chosen by the system will be addressed in the next section. This 
will serve as a useful illustration, in a simple setting, of a phenomenon which actually 
occurs in realistic models with a large number of Yukawa couplings, such as R-parity 
violating models to be discussed in section 4. 


3 Multiple IRQFPs and Landau Pole free domains 

Let us hrst state the two main points we will be discussing in the next two subsections: 


(I) The IR Quasi Fixed Points discussed in the previous section, and studied in the 
literature, are only isolated points on the boundary of the Landau Pole free domains 
(LPfd) which will be dehned by Eqs.(|^) (or by Eqs. QT^) ). Thus the “rectangle 
domains” given by the Quasi Fixed Points encode part but not all of the constraints 
on the allowed low scale values of the Yukawa couplings. 


(11) In extended Yukawa sectors and depending on the ratios dkj/djj there can appear 
simultaneously several IR Quasi Fixed Points for the very same conhguration of in¬ 
creasingly large initial conditions! Actually, some of these points are repulsive and 
others attractive, depending on the way one approaches the given conhguration of 
initial conditions. 


The relevance of point (11) will become manifest through the specihc examples of the next 
sections. We just anticipate here two important and mutually related consequences: (i) 
the attractive IRQFP’s lead to a dynamical suppression of some of the Yukawa couplings 
and (a) the perturbativity bound becomes, numerically, substantially different from the 
Landau Pole bound. Furthermore, it should be clear that the multiple solution scenario 
we will be discussing is not to be confused with the fact that one obtains different quasi 
hxed points for different sets of hnite/inhnite Yukawa initial conditions (like for the typical 
examples of “small tan /3” or “large tan /3” hxed point scenarios EE)- 

3.1 LPfd beyond the rectangular approximation 

Hereafter we determine the equations dehning the Landau Pole free domains (LPfd) be¬ 
yond the crude “rectangular approximation”. It is convenient for the discussion to write 
Eq.(l) at two distinct scales t and f°(< t) (thus t corresponds to an energy scale lower 
than that of t^). Eliminating Yfc(O) in the ensuing equations, one obtains the general cor¬ 
relations between the Yukawa couplings at two arbitrary scales, expressed in a bottom-up 
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approach where the Yukawa couplings at the high scale are cast in terms of their initial 
conditions at the low scale t , 


w(t°) 

_ Yk{t)uk{t°,t) 

1 - akkYkit) //o Mfc(P, t)dt' 

(14) 

where we dehne 



Uk{t\t) = 

Ek{t',t) 

(15) 

j.Pfc (1 - /// Uj{dk 

Ek{t',t) = 

Ek{t') 

Ekit) 

(16) 

Since all the Cki coefficients are 

positive, 12 , a system of non trivial constraints 

on the 


values of the Y^s at the scale t follows, 


Yk{t) < 


1 

JtO '^kit , t'jdit 


(17) 


As can be seen from Eq.(^) each of the above inequalities (for k = 1, ...n) involves all 
n Yukawa couplings at the same scale t. The conjunction of all inequalities dehnes the 
allowed region in the n-dimensional parameter space of Yk{t). We should stress here that 
Eqs.fpTl) describe two different, but complementary, constraints: 

- The strict inequalities preclude the Yk from becoming inhnite anywhere between the 
low scale t and the high scale corresponding to Equations ([T7|) are thus necessary 
conditions guaranteeing the perturbativity of the running Yukawa’s between t and 
But they are obviously by no means sufficient to guarantee perturbativity; when going 
closer and closer to the boundary of the domain dehned by Eqs.(^), higher order loop 
contributions to the running should be included already at intermediate scales between 
t and t^, and a deeply non-perturbative regime is expected to set in in the vicinity of 

Thus the LPfd conditions have a physical meaning at the one-loop level only when 
considered as necessary perturbativity constraints. 

- The domain dehned by Eqs.(^) has yet another meaning beyond the LPfd conditions. 
This is simply the positivity of the W which are by dehnition the squares of the Yukawa 
couplings appearing in the Lagrangian^. (The reader is referred to Appendix A for tech¬ 
nical comments about the derivation of Eqs.(p!7D.) 

Without the positivity constraint, but requiring a perturbativity constraint |Yfc(t°)| < 6, 

where 6 is sufficiently small, one can in principle allow the domain dehned by the reverse 

^The latter couplings are taken real throughout this paper. Note also that for the sake of generality 
we avoid to consider Eqs.(p^ as a requirement for the reality of the Uk functions since such conditions 
are affected by the model-dependent values of the various ratios akj/ajj in Eq.(p^. 










inequality in Eqs.([T7|). With the positivity constraint and YkiiP) < 5 one obtains the 
improved bounds 


Yk{t) < 


(18) 


Ukif, t) +5 akk /to Uk{T, t)dT 

which lead back to Eqs.(p!7D in the limit of large 5 . It is customary to take a/Itt as the 
perturbative bound on the Yukawa couplings hk at the GUT scale, that is 5 ~ l/(47r) when 
normalized as in Eq.(|^). However, one can still get back approximately Eqs.(pTD when akk 
is sufficiently large. This is typically what happens for the top quark Yukawa coupling in 
the standard model or the MSSM: the Landau pole bound and the perturbativity bound 
come up numerically very close to each other even though they have theoretically different 
meanings. In the present paper we exhibit cases where this is no more true. 

Let us illustrate the issue with the simple model of section 2.3. In Fig.l (a) we sketch 
the case |a|/|6| < 1. Point B denotes the unique IRQFP given by the hrst equation in 
(|^) at a given low scale t in the parameter space (Y2(t), Yi{t)). The two thick-thin dashed 
curves correspond each to one of the boundaries delimited by one of the two equations 


Yi{t) < 

Ut) < 


1 

On Ito Ui{t\ t)dt' 
1 

022 //o U2{t',t)dt' 


(19) 

( 20 ) 


together with ([T5| , ^) wherein k = 1,2 and (respectively) j = 2,1. The conjunction of 
Eqs.(pT|, 1^) gives the LPfd domain dehned by the thick curves. (As they stand, the above 


equations are highly implicit in the variables (Yi(f), Y 2 (^))- Specihc forms within some 
approximations have been given elsewhere 


A point lying on one of the thick dashed lines corresponds to a conhguration where 
only one of the two inequalities (^, is saturated, that is only one of the two initial 
conditions (Yi(f°), Y 2 (^*^)) is inhnitely large. It is point B, the unique intersection of the 
thick dashed lines which corresponds to both initial conditions becoming inhnite, i.e. to 
the I.R. Quasi Fixed Point. This illustrates our point (I) in a specihc context. The shaded 
rectangle would be the allowed region if one relies only on the I.R. Quasi Fixed Point. 
The point here is that one is still allowed to go above the IRQFP value for one of the 
Yukawa couplings provided that the other Yukawa coupling remains sufficiently smaller. 
Thus a complete analysis of the theoretically allowed domain and its comparison with 
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Figure 1: The IRQFP configurations for the two-Yukawa system defined in section 2.3. 
The thin lines represent the flow of (Yi, Y 2 ) at the low scale t, when the initial conditions 
(Yi{t^),Y 2 (t^)) are increasingly large. The different flows correspond to a different fixed 
ratio between these initial conditions. 


the eventual experimental bounds should take into account the correlation between the 
Yukawa couplings. Finally we stress that the thick dashed lines are never reached, even 
in the limit of infinite since the only limit is point B. (see appendix B.2 for 

further discussions). In Fig. 1 (a) we show typical trajectories of the Yukawa couplings 
at the scale t. When the initial conditions at t^ are taken large (with any fixed finite ratio 
Yi{t^)/Y 2 {t^)), all the trajectories are attracted towards B, even though they can still go 
outside of the rectangle domain. Although only qualitative and restrained to the case 
of two Yukawa couplings, the above features are quite generic to an arbitrary of Yukawa 
couplings, provided that Eq.(|l^) has a unique solution. This is the case for the MSSM 
(including all Yukawa couplings) |T^ and for the (M+1)SSM The next section is 

devoted to the intricate configuration when Eq. (|T^) allows for more than one solution. 


3.2 Multiple Landau Poles 

We come now to point (11), first noticed in |^. Here we give a general account, illus¬ 
trate the phenomenon in the context of the simple toy model of section 2.3, and leave 
to section 4 the discussion of the more complicated but physically interesting case of the 
MSSM with R-parity violation. As was pointed out in section 2.2, the equation which 
controls the asymptotic low energy behaviour of the running Yukawas when their initial 
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conditions are large, is in general not guaranteed to have only one solution. To start with, 
such a multi-solution conhguration could seem mathematically inconsistent and contra¬ 
dicting the necessary unicity of the solutions of Eq.(2)^ This is actually not the case: 
the multi-solutions are due to the strict inhnity of the initial conditions, and the unicity 
of the running is in fact restored for large (but hnite) such conditions. However, the nov¬ 
elty is that the multi-solutions will all play a role here, as they will correspond either to 
attractive or to repulsive Quasi Fixed points. Actually, contrary to the case of a unique 
solution where the (hnite) ratios = Yj{t^)/Yk{t^) could be varied at will with no effect 
on the quasi hxed point, here the values of Vjk trigger the quasi hxed point towards which 
the system will evolve. 


Let us illustrate this phenomenon with the system dehned by Eq.(^). As was shown in 
section 2.3 and Appendix B.l, there is potentially three different I.R. quasi hxed points, 
i.e. three solutions to Eq.(|^). Only one of these solutions corresponds to both asymptotic 
powers pi,P 2 being smaller than one, while the two other solutions have respectively 
Pi < 1,P2 > 1 and Pi > l,p 2 < 1. From the discussion of section 2.2 one then expects a 
quasi hxed point where both and are non vanishing, and two other quasi hxed 
points where one (and only one) of them vanishes, see Eq. (0. The question now is which 
of these quasi hxed points will be ehectively operating at low energies? We will show that 
the quasi hxed points with one vanishing are the ehective (attractive) ones, while 

the quasi hxed point with the two non-zero is repulsive. This is illustrated in 

Fig.l (b), where now the two boundary lines dehned by Eqs.(^, |^) intersect in three 
distinct points corresponding to the three solutions, and the LPfd is delimited by the thick 
dashed lines. To understand the pattern of the how depicted in the hgure, let us consider 
a point lying close to the boundary of LPfd. Note hrst that, similarly to the case where 
|a|/|6| < 1, this point cannot reach the LPfd boundary (when the two initial conditions 
Yi(t^),Y 2 (t^) become inhnite) elsewhere than on points A, B or C. (See appendix B.2 
for a more technical discussion.) In order to guess the dynamical behaviour close to the 
boundary it is useful to re-write Eqs.(|^, ^) in the form 


Yk{t) = 




(^kk ftO '^k{t , t'jdit 


{k = l,2) 


( 21 ) 


where > 0. This allows to write Eq.(|T^ in the form 

"^Indeed, the running range for t does not cover any gauge coupling Landau Poles so that the system 
defined by Eq. (® still satisfies locally a Lipschitz condition [^, whence the unicity of the solution. 
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Uk{t',t) = 


Ek{t',t) .n, {jloUj{T,t)dTY>‘E^n 

jPfc (/iO Uj{T, t)dr + ej //, Uj{T, t)dTY'^il^n 


( 22 ) 


The important point here is that the denominator in the above equation^ can potentially 
become very small, thus Uk{t',t) very large, when ej <C 1 and simultaneously t' ~ t^. 
But this is exactly the critical region we need to study: the behaviour of Uk{t',t) for 
would determine, through Eq. f^ll) , whether Yk{t) goes to zero when —>• 0, that is 

whether the system is attracted towards points A or C, rather than B, when approaching 


the boundary of LPfd. Because in Eq.(^^ the u^’s are dehned iteratively in terms of 
one another , we are lead to considerations similar to those of section 2.2 where now in 
Eq.(D is replaced by e^. A more sophisticated treatment leading to equations similar 
to Eq■ (pTl) is required though and is deferred to appendix B.2. The bottom line is as 
follows: if we consider the hierarchical regime such as 62 ei < 1 which corresponds 
to a point lying very close to the LPfd boundary between points B and C, then the 
solutions to Eqs. (p.2.7| , |B.2.8|) imply unambiguously that //o ui{t\ t)dt' ~ 62 where | A| 
is a positive number uniquely determined by these equations. Combined with Eq.(^), 
this behaviour means that the system is dynamically attracted to point C. Similarly, 
if we consider ei <C 62 < 1, then the system will be attracted to point A. This leads 
altogether to the flow pattern of Fig.l (b). Actually, it is still possible to land on point B, 
but only at the price of a delicate hne-tuning corresponding to the trajectory satisfying 
Yi{t) = especially when |a|/|&| is much larger than one. In Fig.2 we show a purely 

numerical analysis staring directly from the differential equations ( 0 ) with a = l(resp. 6 ) 
and b = 6 (resp. 1). The results conhrm nicely the above theoretical arguments. Fig.2 (a) 
illustrates the insensitivity to the values of ri 2 = Y^/Y 2 near the IRQFP. In contrast, 
one notes the high sensitivity to this ratio in the conhguration of Fig. 2 (b), due to the 
presence of multiple Landau Poles. In particular, the attraction to point B in Fig. 2 (b) 
occurs at the price of hne-tuning ri 2 to 1 to less than 0.1 per mil, otherwise the system is 
strongly attracted to points A 01 C depending on whether ri 2 is greater or smaller than 
one. 

To summarize this section on the two-Yukawa system, we have shown in general that when 
ri 2 is neither inhnite nor zero, the LPfd boundary closes only on a set of isolated points, 
point B in the case of Fig.l (a) and points A, R, C in the case of Fig.l (b). Most of the 
LPfd boundary, the thick dashed lines, can be closely approached but is actually never 
reached, the system being eventually driven to one of the above points. In Fig.l (a), only 

^ In all models we consider, cikj/'^jj is a positive number; for the two-Yukawa case of section 2.3 we 
took akj/ajj = |a|/|6|. 
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1^2 1^2 




Figure 2: A numerical illustration of the behaviour of the toy model Eqs.dl^) with |a| = 
1, |6| = 6 in a) and |a| = 6, |6| = 1 in b). The running Yi{t),Y 2 {t) are evaluated at 
t = 66. The smallest initial condition Yf is varied from 0 to 1000. Near the IRQFP’s, 
the insensitivity to ri 2 in a) and the high sensitivity to ri 2 in b) are manifest. 


when ri 2 is infinite or zero (i.e. one of the initial conditions remains finite) does one retrieve 
the effective fixed point a la Hill, in which case point B becomes located on one of the 
two axes. Of particular importance is the case of Fig.l (b). It shows how the existence of 
multiple IRQFP’s fits nicely with the phenomenon of dynamical attraction to or repulsion 
from such points, leading to the unusual suppression of some couplings. This is quite a 
general mechanism which occurs beyond the two-Yukawa system as will be illustrated in 
the next section. Finally we would like to insist on the fact that although a perturbative 
analysis is trustful only far from these IRQFP’s, the latter remain nevertheless very useful 
in understanding the generic behaviour of the low energy perturbative couplings, as they 
play the role of a far away beacon for such a behaviour. 


4 ^„-MSSM 


In this section we concentrate on a phenomenological application in the context of the 
MSSM with R-parity violation. The unusual behaviour near Landau Poles was indeed 


first noticed in this context 14 


Let us first recall the form of the superpotential in .^p-MSSM: 


13 









Wf,^ = W + W^ + 

where W is the R-parity conserving part 


(23) 


= (hEhLlH.Bjj + {hohQl.H.D^E + (l‘uhQ'L-H2Ui + IJ.H 1 .H 2 (24) 
Wji induces lepton number violation, 

Wf = + tJ,Li.H2 (25) 

and Wp induces baryon number violation, 

= ^K^iU'E&ED'k. (26) 

The superhelds Ll,Ql^Er,Dr^Ur denote lepton and quark SU{2) doublets and anti¬ 
singlets, Hi,H 2 the two Higgs doublets. Summation on SU{3) color indices is implicit 
and the dots {A - B = eabAaBb) dehne SU(2) invariants The i,j, k = 1, 2, 3 are generation 
indices, and summation is understood for repeated indices. 

We will rely hereafter on ([^) to write down the RGB’s at the one-loop level. The 
relation between the above A’s and their matrices A is as follows: 


\jk — i.^E^)iji ^ijk — {■^D>‘)iji ^ijk — {■^U^)jk- 


(Note that the notation is consistent with that of |^|, but differs for instance from |R 
only in A".) The matrices A^k and A^k being anti-symmetric, lead to a total of 18 
couplings, while Aj^k give 27 extra couplings. We will not consider hereafter the bi-linear 
couplings K* and /x as their renormalization group evolution does not affect that of the 
Yukawa couplings we are interested in, 0 , m- 

The issue of theoretical bounds on RPV couplings has been addressed in various papers 


151 , hq , pTl, |27|, ^ , either from the point of view of exact infrared fixed points, when they 


exist, or from that of perturbativity and/or effective hxed points related to Landau Pole 
considerations. However, often a rather reduced number of couplings has been considered 
which may not exhibit the behaviour we described in the previous sections. Typically, 
the larger the number of couplings, and accordingly the dimension of A4 in Eg. ([I0|) , the 
more likely it is to find asymptotic powers significantly greater than 1 as solutions to 
this equation, leading to a dynamical suppression of some Yukawa couplings due to the 
Landau pole. 


As was first shown in |jT^, the system of 6 couplings ht, hr, A 233 , A 333 , A 323 does ex¬ 


hibit such a behaviour where A 333 decreases while the other couplings increase towards 
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their perturbativity bound. 


On the other hand, there are also numerical studies of the simultaneous renormalization 
group evolution of all 45 Yukawa couplings, putting the accent on the correlation between 
experimental bounds at low energy and initial values at the grand unihed scale pS . 


Here we take up a fairly complementary approach. We consider various sectors with 13 
couplings chosen in such a way as to encompass previous studies as subclasses and in the 
same time to be realistic enough for phenomenology. These sectors allow also a direct ap¬ 
plication of the analytic criteria developed in the previous sections and a clear numerical 
illustration of these criteria. 


RPVi3 sectors: Each of these sectors has 13 couplings and is dehned as follows. For hxed 
k and hxed l{^ k) we consider the RGB’s associated to the system 

RPVi3 = {ht, hb, hr, Xkiq, Kkk, Kki, Kik, Kkk, Kki)^ wliere g = 1, 2, 3 (27) 

Thus one has for each case three A’s, three (A")’s and four (A')’s. The RGB’s in the RPV 13 
sectors are explicitly written in appendix G. One can see that they have strictly the form 
of Eq.(H) when k ^ 3. For k = 3, Eqs.(G.2, G.3, G. 6 , G.IO) deviate by one extra term 
from this form, but as we will see throughout the numerical analysis, the behaviour is still 
dictated by our general criteria. Actually the specihc choice of the set (^TD is motivated 
by the opportunity to remain as close as possible to the form of Eqs.(2). We stress here 
that this is not necessarily the only possibility, but we stick to it as a working hypotheses 
which allows for general enough, yet analytically tractable, illustrations. 

4.1 Asymptotic powers 

To start with, it is instructive to follow the variation of the asymptotic powers when the 
number of active RPV couplings is increased. 

RPVg: Gonsider the set of non zero couplings {ht, hb, hr, X233, X'^33, X'^23)- This is a 
subset of Eq.(p7[) with k = q = 3 and I = 2. Here the order of the antisymmetric indices 
is irrelevant, since the evolution equations in the RPVe subset are blind to the sign of the 
couplings, as can been seen from Appendix C. In this case Eq. ®) reads 
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( 28 ) 


Solving this equation in the regime where all 6 's = 1, one hnds that Pt,Pb,P233,P323 < 1 
while P 333 = I and p,- = || meaning that the running A 333 and hr at any sufficiently low 
energy scale will decrease to zero with increasing initial values at the GUT scale. However, 
since p'^^^^Pr are close to 1 this decrease is relatively slow, [Q. Actually Eq. p 8 |) has yet 
another solution with all asymptotic powers smaller than 1 ; but this solution turns out 
to be repulsive in the sense of sections 3.2 and B.2. In Fig. 3 we illustrate the evolution 
of the RPVg system, solving purely numerically the corresponding RGB’s as given in 
appendix C. 


A* 



Figure 3: The running couplings of the RPVq system evaluated at t = 66 ~ 

Log[MQjjrp/Mz], as a function of their common initial condition A°. A 333 and Xr decrease 
very slowly even for A° S> Xp^^f 

The decrease of A 333 and hr with an increasing common value of the initial condition 
A° nicely conhrms the above expectations. However, since the corresponding asymptotic 
powers are close to 1 , the dynamical suppression is too slow to be sizeable before the 
perturbativity bound at the GUT scale, dehned for instance as A^g^^ = ~ 3.54, is 

reached. 
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RPVg: Adding A 233 , A 323 , A 332 to the previous set of couplings and solving the corre¬ 
sponding Eq.([T^) with the help of algebraic packages (the relevant 9x9 matrix Afg 
can be straightforwardly extracted from Adis in Eq.(p9D below) one now Ends a larger 
value for 1.67). The decrease in A 333 becomes strong enough to build in signih- 

cantly before the perturbativity bound at the GUT scale is reached. Note that there are, 
here too, more than one set of solutions, actually three, but all of them have the above 
value for P 333 . The dynamically chosen solution has also pb — 1.16 and P 323 — 1.3. Again, 
we have illustrated in Figs.4 (a), (b) the numerical results obtained directly from the 
RGE’s of the RPVg system and which conhrm the above asymptotic power behaviour. 
Gomparing Fig.3 with Fig.4 (a), we see how the activation of a bigger number of RPV 
couplings strengthens the bounds on some of these couplings, in particular the drastic 
drop of A 333 in the RPVg as compared to the RPVg case, or to the case of smaller number 
of active couplings like in |^. The price to pay is a looser bound on the extra couplings 
as shown in Fig.4 (b). 

A?; A,; 




Figure 4: The running couplings of the RPVg system evaluated at t = 66 ~ 

Log[MQjjrp/M"^], as a function of their common initial condition A*’. 

RPV13: Addition of new couplings does not necessarily make the dynamical suppres¬ 
sion of one single coupling more pronounced, but can rather propagate this suppression 
to other couplings. This is indeed the case when all 13 couplings as dehned in Eq.fpT]) 
with k = 3, are simultaneously activated. The corresponding matrix of Eq. (pd]) is easily 
read from appendix G: 
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(29) 


In this case (again with all 5 's = 1 ), one hnds, with the help of algebraic packages p 5 
14 different solutions to Eg. (p! 0 D . All these solutions have in common P333 > 1.43 and 
P332 ~ In particular, the dynamically chosen solution exhibits new features: P333 has 
decreased slightly to ~ 1.46 as compared to the RPVg case, but now P323 and P332 are 
1.3 and P233 1 . 06 ; the asymptotic powers of all other couplings remain smaller than 

one. This means that the three couplings A323, A333, A332 are dynamically driven to zero 
by the Landau Pole, and in any case decrease significantly, well before the perturbativity 
bound is reached, as confirmed by the numerical analysis shown in Fig .5 (a), (b); A233 is 
also dynamically driven to zero but decreases very slowly since its asymptotic power is 
very close to 1. 

[Note that all the above conclusions hold as well for I = 1, i.e. for A133, A333, A333, etc... 
since Afi3 is / independent. ] 


4.2 Theoretical bounds on the RPV couplings 

As can be seen from Figs.3,4,5, the couplings which increase with A° actually saturate 
and become quickly weakly sensitive to A°. This is why the perturbativity bound for the 
low energy couplings is sometimes associated to the Landau Pole. This is indeed typically 
the case for the top, bottom and r Yukawa couplings in the R-parity conserving MSSM. 
However, the fact that the perturbativity bound and the Landau Pole are in general dis¬ 
connected is nicely illustrated by those couplings which decrease with increasing A° and 
are sensitive to its values. For instance, in RPVg and RPV 13 , A 333 drops by more than 
a factor of 2 of its maximal value, when A° reaches a/ 4^- Moreover, the perturbativity 
bound on A 333 (taken here at Mg^sy = 1 TeV) drops drastically, from 0.52 in RPVg to 
0.12 — 0.13 in RPVg and RPV 13 , as can be seen from the figures. These bounds are way 
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Figure 5: The running couplings of the RPVis system evaluated at t = 66 ~ 

Log[MQUj,/M‘^], as a function of their common initial condition A°. 


below the ones usually quoted. 


In order to have a realistic comparison, also with the existing experimental bounds, one 
has to require consistency with the experimental top, bottom and r fermion masses. This 
translates into the following constraints 


mt(2y2GF)2 mb,r{mt) {f2^/2GFY^ /Qn^ 

= - —5 -, hb,r{mt) = - 3 - (30) 

sm p cos p 

where Gp is the Fermi constant (= 1.166 x 10“^ GeV), mt is the running top quark mass 
at its hxed point value mtijnt) = mt, and mb^rimt) are the running bottom and r masses 
at this value. We take mt = 165 GeV, mbimt) = 2.76 GeV and mrimt) = 1.78 GeV. Fur¬ 
thermore, we will also require unihcation of the three gauge couplings at Mqut = 2 x 10 ^® 
GeV. Keeping in mind that we are not seeking here a very precise gauge unihcation anal¬ 
ysis since the running is performed only at one-loop level, we tolerate a deviation of the 
value of as{Mz) from the one determined by LEP precision measurements, as well as a 
susy scale ~ Mz in order to achieve one-loop gauge unihcation with the proper normal¬ 
ization ((5/3) X gf{GUT) = g^lGUT) = g‘^{GUT)). It should be clear, however, that a 
more sophisticated treatment of this sector would not change the main features related 
to the dynamical behaviour of the RPV couplings. 
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We take the following values for the hne structure constant Oem, and sin^ 9w at the 
Mz scale 


aem{Mz) = 1/127.938, = 0.1169, sin"^ 9w{Mz) = 0.23117 


(31) 


(of which as deviates by 0 . 2 % from the experimental value 


In order to study the effect of Landau Poles, the numerical algorithm will have to meet 
the requirement of maximizing the initial Yukawa conditions at the GUT scale, while 
keeping consistency with Eq. (|30|) . This is not straightforward due to the large number of 
couplings. As we stressed in previous sections, Landau Pole constraints (or for that matter, 
perturbativity constraints) delimit hypervolume domains in the space of couplings. The 
numerical values obtained correspond only to specihc directions in this space, depending 
on which set of couplings is made large at the GUT scale. It is important to determine 
the directions which allow the largest possible values for all the couplings. 

In practice, we proceed as follows: for a given tan/3, and assuming a common susy 
scale much higher than the top mass, we hrst run the (non supersymmetric) standard 
model top, bottom and r Yukawa couplings from their values at m* scale as given by 
Eq.(^) up to Msusy scale. From there we run these couplings up to the GUT scale 
within the (R-parity conserving) MSSM. There we switch on all the RPV couplings of a 
given RPV„ sector, giving them in a hrst step a common and large value at the GUT 
scale ycuTi ^^e n Yukawa couplings of the RPV„ sector down to Msusy, re-adjust 
ht, hi,, hr to their previously determined values consistent with Eq.(^ 0 D and run the RPVn 
sector back to Mgut- If at least one coupling is greater than the perturbative bound 
(taken to be then the procedure is iterated starting from a smaller ia the 

opposite case Xqut is increased and the procedure iterated until the largest coupling at 
the GUT scale reaches the limit (within 1 per mil). However, even then, there is no 
guarantee that the other couplings have been optimally maximized, since the procedure 
starts with a unihed value of the RPV’s at the GUT scale. This is why a further step 
is taken after numerical convergence: start with the obtained solution at Msusy, scale up 
some of the RPV’s and down some others, by the same factor, then run up to Mgut 
and check for perturbativity. Again the procedure stops after convergence to the optimal 
scaling factor. The numbers thus obtained at the susy scale correspond to theoretical 
perturbativity bounds on the various couplings. It is important to note here the existence 
of a seesaw effect between some specihc sub-classes of RPV couplings. The above inverted 
scaling works when these sub-classes are properly chosen. For instance, in RPV 13 , the 
dynamically suppressed couplings X 313 , A 333 , A 33 ; (see Fig.5) can be increased at the price 
of lowering the remaining ones. 
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In Fig .6 we illustrate these effects as a function of tan (3 in the case RPV13. The difference 
between Figs, (a), (c) on one hand and Figs, (b), (d) on the other, shows clearly that 
it is possible to increase, for instance, the magnitudes of As^a, A 333 for any value of tan/3, 
but various other A', couplings should be decreased accordingly, in order to remain 
consistent with the constraints described above. This seesaw effect is actually far more 
general than what is illustrated in the figure and occurs for the sets of A’s, (A')’s and 
(A")’s either separately or in some combinations of them; it is triggered by the fact that 
not all couplings can be simultaneously increased (or decreased) without conflicting with 
the dynamical suppression mechanism explained before. Also the magnitude of the effect 
tends to increase with decreasing tan /3. For instance, if tan /3 < 10 one can raise A 33 ; up to 
~ 0.97 at the price of reducing all the remaining RPV13 couplings down to < 0.17 for the 
A’s and A 333 , and down to < 0.3 for the other A', A". For higher tan j3 values these bounds 
drop further down: when tan/3 ~ 60, A 33 ; ~ 0.84, while A 33 ;, X 231 ~ 0.15, X 313 ~ 3 x 10“^, 
X 311 , X 312 , AJ33 ~ 7.7 X 10“^, A333 ~ 6.8 X 10“^ and Ag^g, A33; ~ 0.12. (Similar effects obtain 
if one chooses to raise the bound for A 333 ; we will come back to that in the next section 
when comparing with the present experimental limits.) We thus see that Fig .6 is with 
many respects a conservative illustration of the perturbativity bounds one can obtain. 

Note that although the top/bottom/r Yukawa couplings are uniquely determined in 
terms of tan /3 through the fermion masses at the relevant low scale, Eq.(^), their running 
values above the susy scale are still affected by the RPV system. In particular, the GUT 
scale values of hfe, hr remain small for tan /3 < 40 where A 333 or X 313 are those which satu¬ 
rate the perturbativity bound, while for 40 < tan [3 < 60, hb then hr take over gradually 
in saturating this bound. 


The main lesson to draw from this section is the importance of including as many couplings 
as possible in order to achieve a complete study of the perturbativity bounds. The results 
obtained here for the RPV13 system yield much stronger perturbativity bounds on X' and 
A" than previously found ||2^, (see for instance the bracketed numbers displayed in the 


compilation Table 1 of ref. p^). The latter vary in the range 1.04 — 1.23 while ours can 
hardly reach 1 for some couplings while the others are anyway bound to be much smaller. 
The fact that stronger bounds ensue from activating a larger number of couplings was 
already exemplified in [^, |^, since there the upper bound 1.23 was on \fj2i ^ 331 ^ rather 
than on individual couplings. A similar effect is actually observed in our numerical study, 
nonetheless it should be considered as a “phase space” effect which is complementary, but 
distinct, from the dynamical suppression effect we have studied throughout the previous 
sections. 
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4.3 Comparison with the experimental bounds 


Bounds on RPV couplings from experimental data have been extensively studied in the 
literature, ranging from nuclear and atomic physics to astrophysics, as well as high energy 
accelerator physics, (see [Q, and references therein and there out). Since we are 


considering several A, A', A" in a time, a systematic comparison with the above bounds 
is certainly out of the scope of the present paper. Indeed, these bounds are usually 
established keeping a minimal set of couplings in a time |] 30 | , | 3 l|, in particular for the 

very strong ones on the products A x A" from proton stability or on A x A, A' x A' 

from neutrino masses . We shall thus limit the discussion here to some representative 
examples. We stress here that the various experimental limits on individual couplings at 
the weak scale (see for instance Table 1 of ref. [^), if one assumes a common susy mass of 
rfi = 1 TeV, can be in general easily overwhelmed by the theoretical perturbativity bounds 
we have found in the RPV13 system. Comparing for instance to the conhguration of Fig.6, 
the experimental limits on A321, A322, A323 ~ 0.7 are about a factor of two weaker than our 
perturbativity bound, except for A323 at low tan /3 in the conhguration of Fig.6 (b) where 
they become comparable. However, in the latter conhguration the perturbativity bounds 
on A233, A323, A332 < 0.3 are stronger than the experimental limit^. Of course there are 
conhgurations where the experimental limits become, at least for some couplings, stronger 
than the perturbativity bounds. This is true in particular if one takes into account recent 
limits from proton stability [Q and neutrino physics [^. If we pick up for instance the 
conhguration where A33; ~ 0.84 (see previous sub-section) then the perturbativity bound 
on the products X313X331, X3iiX'^:^i, X3i2X'^;^i ~ 2.5 x 10 “^ — 6.5 x 10 “^ is overwhelmed by the 
proton stability limits by several orders of magnitude 1 ^^. For the same conhguration one 
has the perturbativity bounds A333 ~ 6.8 ( 3 . 6 ) x 10 “^ and X313 ~ 3 ( 1 . 9 ) x 10 “^, where the 
numbers in brackets correspond to a lower susy scale mgusy = 100 GeV. These bounds are 
clearly stronger than the experimental limits obtained from specihc r lepton and Z boson 
branching ratios measurements, but they remain two orders of magnitude weaker than 
those from neutrino physics [Q. However they are still meaningful in various respects: 
(i) they correspond to a situation where many couplings are active simultaneously, while 
the bounds from neutrino physics and proton stability are obtained within reduced sets 
or simplifying working assumptions ^, 3 ^, barring possible cancellations ehects which 
could in principle weaken substantially the latter bounds, (ii) They are of the same order 
both for rrtsusy = 100 GeV or 1 TeV, while experimental bounds from neutrino physics are 


®Note that even though we have evaluated the running couplings at 1 TeV, these values are supposed 
to be frozen and remain valid at the lower scale due to the simplifying assumption in this particular 


illustration that the susy spectrum is at or above nig 


= 1 TeV. 
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derived assuming rususy = 100 GeV. (ui) In both proton stability and neutrino physics it 
is not clear how low energy QCD effects between virtual quarks and squarks at and 
vertices would affect the extraction of the experimental bounds. The running RPV’s are 
in principle free from such uncertainties insofar as their low energy values are frozen at 
the susy scale. 

In the case where mgusy — 100 GeV, most of the 2a limits on the RPV’s are a factor 
of ten stronger than for rUsusy — 1 TeV. The perturbativity bounds are then overrun for 
small and moderate tan/? values, in particular for the A’s of RPV13, but not necessarily 
for the (A')’s and (A")’s. For instance, tan/3 ~ 10 leads, for a conhguration where is 
maximized (~ 1.09), to A^'g; ~ 0.34, A'^, ~ 0.28, A'333 ~ 0.18, - 0-3, while 

~ 0.18. For tan/? as large as 60, A'/g^ maximizes to 0.62 and A23; — 0.15, A33; ~ 0.11 
•^333 — 4 X 10 Agg^, Ag^g — 0.11, AJgg ~ 5 X 10 Ag/i, A3;2 — 7 X 10 ^ and Agjg ~ 2 x 10 
There is a tremendous drop in all couplings which makes the above bounds stronger than 
most of the experimental limits (compare Table 1 ref. [^] with rUsusy ~ 100 GeV). They 
remain however much weaker than the limits from neutrino physics and proton stability, 
even though some of them (e.x. A322) come very close to these strong bounds, within a 
factor of two. Taking these limits at face value (albeit points (i) - (in) made above), they 
thus favour scenarios with large tan/?, which would anyway be what is typically needed 
in order to have part of the susy spectrum in the electroweak scale range. 

We should keep in mind, however, that the above discussion has been carried out within 
specihc conhgurations of the initial GUT scale conditions for the various RPV’s, where 
a common value was assumed, at least as a starting point for the numerical algorithm. 
Ghoosing hierarchical initial conditions would allow to suppress a larger number of cou¬ 
plings at the electroweak scale, as can be seen from the general structure of the solutions 
Eqs.(|, I), and thus to meet eventually with stronger experimental limits. 

To summarize, we have shown in what sense providing one single number for the 
perturbativity bound as is sometimes done, is rather restrictive and cannot be the end 
of the story. Indeed, on one hand such a bound can vary sizeably with the number of 
active RPV couplings, and on the other hand, due to the fact that there are perturbativity 
regions rather than perturbativity bounds, much like the illustration of Fig. 1, a seesaw 
mechanism, involving various RPV couplings can be operating due to the presence of 
attractive and repulsive effective hxed points. This hints to the necessity of including 
simultaneously as many RPV couplings as possible both when determining experimental 
limits and when studying theoretically the correlations induced by renormalization group 
evolutions. 
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5 Conclusion 


In this paper we have considered the running of an arbitrary number of Yukawa type 
couplings to one-loop order, in models with an extended sector of such couplings and a 
gauge sector. We studied in particular the general features of Landau Pole free domains 
for the Yukawa couplings and the related Infrared Quasi Fixed Points. We pinpointed 
the existence of new structures which can be interpreted in terms of multiple repulsive 
and attractive quasi hxed points and developed the analytical formalism which allows to 
determine such structures. An interesting consequence is the dynamical suppression to 
zero of some components of the quasi hxed points, a fact which contrasts with the usually 
expected behaviour. 

We then showed that these new conhgurations appear naturally in the context of the 
MSSM with R-parity violation. In particular, an increasing number of R-parity violating 
couplings induces a suppression of some of these couplings when the initial conditions are 
large. A notable example is the suppression of Aua, i^'kkk) ^ski to zero at the Landau Poles. 
This would be an interesting theoretical justihcation of the smallness of such couplings, 
which goes along with stringent experimental limits as well as the necessity of stabilizing 
the proton and the Neutralino LSP, were it not for the fact that in practice one should 
keep far from the Landau Poles for a consistent perturbative treatment. Nonetheless, 
taking into account the perturbativity bounds and constraints from the physical quark 
masses, this suppression mechanism translates into a seesaw effect involving all the cou¬ 
plings. We have demonstrated how this mechanism strengthens the theoretical bounds 
and made a quick comparison with the existing experimental limits. Contrary to what 
one could have naively expected, the simultaneous inclusion of a large number of R-parity 
couplings in the evolution equations, together with the limits from experimental data, can 
lead to more severe bounds on these couplings than in the case where they are studied 
individually or in small sets. This requires that the analyses of experimental limits be 
also carried out while including simultaneously a large set of couplings. 
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Appendix A: Technical comments on the LP condi¬ 
tions 


It is straightforward to see from Eqs. p!5|) that the conditions given in Eq.(p!7[) are 
sufEcient to ensure positivity of the y^’s at all scales between t and The necessity of 
these conditions requires some more care, since there is still the logical possibility that 
some Uk{t',t) and 1 — akkYk{t) It' Ukit",t)dt” may turn negative simultaneously at some 
scale t', < f < t. One could even imagine both those functions taking zero values 

simultaneously thus evading the Landau Pole. We show hereafter that such a situation 
cannot occur. Indeed, starting from Uk(t,t) = 1 (see Eq. (|T^)), if Uk(t',t) —>• 0 for a given 
k and t' between and t, then necessarily |1 — ajjYj{t) //, Uk{T, t)dr| —*• cx) for at least one 
j ^ k as can be seen from Eqs.(|T^ |T^, |). But this can occur only if \ukit,t)\ —>• oo for 
some i in the interval [f,t]. Then in the vicinity of t, the Yukawa coupling behaves 
like 


m 




IiUj{T,t)dT 


(A.l) 


Since ajj > 0, IJT^, and the (large) integral is dominated by the contribution of Uj{i, t), 
then Yj{i) has a negative sign, which contradicts its dehnition. We thus conclude that no 
Uk can vanish in the physical region [t^,t]. 

Finally, to be complete one should mention the logical possibility that the Uk functions 
display a non continuous jump at some scale and change sign without going through zero. 
In this case the positivity of the YkS would require the denominator l—akkYk{t) ft' Uk{T, t)dT 
to change sign too during the jump; since this change remains continuous in t' the de¬ 
nominator goes necessarily through zero. Such a conhguration exhibits a Landau Pole 
only at intermediate scales between and t. Thus it does not £t to the usual physical 
pattern where a Landau Pole appears only at the highest energy scale of the problem, 
signaling the need for new physics in the vicinity of that scale. We thus disregard this 
mathematical conhguration altogether. 


Appendix B 

B.l: Non uniqueness of asymptotic behaviour 

We study the following equation, 

/ Pi \ / a ^ W (1 - pi) 0[1 - pi] 5i \ 

\P2 ) V 7 / V “ ^2) 0[1 - P2] 82 ) 


(B.1.1) 
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Consider the case = ^2 = 1. 

(1) If Pi < 1,P2 < 1 one must have 


Pi = 


P2 = 


a + j3 + a6 — jS'j 
(1 + a)(l + 6) - 13'y 
-y + 6 + a6 — P'y 
(1 + a)(l + (5) — /37 


< 1 


< 1 


(2) If Pi > 1,P2 < 1 one must have 


(B.1.2) 

(B.1.3) 


Pi = 


P2 = 


/9 


1 + 5 

5 

TT5 


> 1 


< 1 


(3) If Pi < 1,P2 > 1 one must have 


(B.1.4) 

(B.1.5) 


Pi = 
P2 = 


a 


1 + 0 
7 

1 + 0 


< 1 


> 1 


The necessary and sufficient conditions for the simultaneous realization of 
and (3) 


(B.1.6) 

(B.1.7) 

cases (2) 


l3>l + 5 , 5>-l (B.1.8) 

7>l + o , o>—1 (B.1.9) 

are found to be necessary and sufficient also for case (1). These equations are thus the 
criteria for the existence of three solutions for Eq. (|B.l.l|) . Note that there are no solutions 
where pi,P2 are simultaneously > 1. 

B.2: Attractive/repulsive IRQFP’s 

Let us now consider the two-Yukawa system of section 2.3. In this case o = 5 = 0 
and /3 = 7 = |a|/|6| in Eq. (|B.l.l|) , so that when |a|/|6| > 1 there will be three different 
I.R. quasi fixed points corresponding to the three configurations discussed in the previous 
subsection, while |a|/|6| < 1 implies a unique solution associated with the configuration of 
case (1). Here we wish to complement, somewhat technically, the qualitative discussion 
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of the two-Yukawa system carried out in sections 3.1 and 3.2. [We skip, though, for the 
sake of simplicity, the discussion of the critical value jlj- = 1.] 

We discuss hrst the nature of the boundary of LPfd illustrated in Figs.l (a), (b). The 
thick dashed lines in these hgures, although representing this boundary are actually never 
reached, even when the two initial conditions are strictly inhnite, except in the points 
A, B or C. Indeed, in any other point on the thick dashed lines, only one of the two 
inequalities in Eqs.(|T^, would (by dehnition) saturate, implying that one and only 
one of the two initial conditions should be hnite. In this case the structure of Eqs.(§, H) 
shows easily that (independently of the magnitudes of the asymptotic powers pi) one of 
the two has to vanish, which contradicts the fact that the point we consider was 

supposed to be distinct from points A or C. 

As was discussed in section 3.2, it is crucial to study the behaviour of Uk{t\ t) as a function 
of t' in the vicinity of t^. [The scale t at which the IRQFP’s are evaluated remains hxed 
and far from Here we give some elements of the derivation of the equations which 
control this behaviour. The ongoing arguments are actually valid for an arbitrary number 
of Yukawa couplings even though presented for simplicity in the context of the two-Yukawa 
model. In view of Eq. (p^ , one expects without loss of generality the following behaviour 
of Uk in the vicinity of 


Uk{t' ~ t°,t) ~ K — (B.2.1) 

^3 

where the ej > 0, and the aj are some real numbers which we wish to determine. For 
simplicity we will assume that there are just two types of e^’s: inhnitesimal ones, which 
we then take all equal and denote by e, and non inhnitesimal ones. Thus, depending on 
the regime we are interested in, taking 6^ = 6 for some fc’s means that we consider a point 
close to the LPfd boundary lines dehned by Eq.(|2l|) for this specihc set of Y^’s, but far 
from the others. 

In order to handle the implicit dependence on e in Eq. (|^) correctly, one has to split 
//o Mj(t, t)dT in the following way: 


/ uATA)dT = 

Jto Jt<- 


P+An /•P+At2 /-P+Ats rt' 

+ + + 

P Jp+Aii Jt0+At2 Jto+Ats 


(B.2.2) 


where 


Ati -C e, At 2 ~ e, e -C Afg -C 1 (B.2.3) 

Indeed, within an expansion in the small parameter e each one of these regions will 
contribute differently to the powers aj. Most generally one should then consider three 
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different sets of aj, i.e. 


Uk{t^ + AtA) ~ ( 1 ) 

if 

At < Ati 

(B.2.4) 


if 

Ati <At< At2 

(B.2.5) 


if 

At2 <At< Atg 

(B.2.6) 




Note that there is no contribution associated to the last term on the right hand side of 
Eq. ( [B.2.2|) since r ^ in the corresponding integration region so that no 1/e behaviour 
is expected for Uk{T,t) there, see Eq. (|^) . Also one can always choose, without loss of 
generality, a Ati sufficiently smaller than e so that the first term on the right hand side of 
Eq.( [B.2.^) can always be neglected and there is no contribution from Eq.( [B.2.^) . Using 
Eqs. (p.2.2HB.2.6|) in Eq.(^) one obtains, after some straightforward but tedious analysis, 
a system of coupled equations for the as follows. 


^( 2 ) = 


^( 3 ) = M-Af 

where the matrix M. is given by Eq. ([Tl|) and A^Ae are column vectors defined by 


(B.2.7) 

(B.2.8) 


^ af 

(4“'*); = ef> + e, (1 - 5,), (x = 2,3) 


(B.2.9) 

(B.2.10) 


with 


0f) ^ 1 - e[af - a®] e[af] af - e\cyf - a®] e[af] af 


(B.2.11) 


0,- = 9[af - af - 1] e[af - 1] (1 - af) - - af + 1] 9[af] af (B.2.12) 

= 1 — 9[af^ — af^ — 1] 9[af'^] — 9[af^ — + 1] 9[af^ + 1] (a® + 1) (B.2.13) 

Here 9 is the Heaviside function, and 5j = 1 (resp. 0), if Cj = e (resp. Cj S> e ). Thus, 
the 5j’s parameterize the various conhgurations of closeness and remoteness from the 
various LPfd boundary hypersurfaces. [For instance, in the two-Yukawa case of section 
3.2, = 0 correspond to the hierarchical configuration e = ei -C 62 -] The 

above formulation is quite general, applicable to any number of Yukawa couplings. We 

















actually tested it with Mathematica, up to the case of 6 couplings, the RPVq system 
of section 4. Here we stick for simplicity to the two-Yukawa system with |a|/|&| > 1. 
Solving explicitly the coupled equations ([B.2.7| , [B.2.8|) in the regime e = ei -C 62 , we 
hud only two consistent sets of solutions, = a® = —{a/hY^af^ = af'’ = |a|/| 6 |) or 
= a® = = 0, = |a|/| 6 |). Taking into account Eqs.( [B.2.2| , [B.2.5| the 

hrst solution implies that in Eq. f^ll) //o U 2 ~ e“l“l/l^l+e^“l“l/l^l+£nite part ~ so that 

Y 2 (^) ~ —>• 0 , and similarly the second solution leads to h 2 (^) ~ 0 , while 

in both cases Yi{t) remains hnite. This proves in general the behaviour shown in Fig. 1(b), 
where the flow is repelled from point B and attracted to point A with a strength increasing 
with |a|/| 6 |(> 1). A similar behaviour occurs for point C in the regime e = 62 ci. 


Appendix C: One-loop RGE’s in RPV13 sectors 


Following we write here explicitly the renormalization group equations to one-loop 
order which govern the evolution of the system of 13 couplings (or less) defined in Eq.(^). 
The following RG equations are valid for one given k and /, A; = 1, 2, 3 and I Y k. It is 
worth noting the extra terms in Eqs.(C.2, C.3, C. 6 , C.IO) when k = 3. These terms violate 
the general form of Eq.(|^) and the insensitivity to the sign of the various couplings. As 
they stand, all couplings are assumed to be positive. A change of sign in one of the 
couplings Xki 3 , hr, hh translates into a sign flip in front of the square roots. 


—-h 3gl —gl - 6hj - hi - {X^kk + ^kki + ^ilk)^k3 - Y-hk3)Xk3k ~ 

16vr^—= -h 3gl —gj - hi - 6hl - h^ - A^^, 2X^ii. + QX^j. 

+2(Aifci + X 2 I 1 + XYl))hk3 - (1 - ^k3){^k3k + 2Afcfc3 2(A3^3 + ^2k3 + '^sD)} 

-^KlkKi3hlhl5k3 


(C.l) 


(C.2) 


16vr^——h^{35f3 + ~ - (A^^^ X\i2 + ^Xh^, + 3X^^i. + dA^^; 3XYk)^k3 

— (1 — <5fc3)(A^33 -|- A^32 -|- 4A^33)} — 3k3 


(C.3) 


^Q'^^-^^kii — ^lii{^9i + ^92 ~ hl{6k3 + his) - 4A 


2 

kll 


4-^M2 


^^kl3 


^i^kkk + ^kkl + ^klk + ^ 


Ikk 


)} 


(C.4) 


,^hYi2 — Xli2{3gl + 3g‘l h^(5fc3 + <^/3) 4A^;]_ 4A^;2 4A^j3 3(A^^^^-I-A^^^^-I-A^f;^-I-A^^^)} 


dt 
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— + 3fi'2 ~ 2(1 + 5k3 + 5i3)hl - AXh^ - A\li2 - 4A 


dt 


2 

kl3 


■3(Afcfcfc + Xkki + ^kik + ^ilk)} ~ ^\J^kk^kis^Thb^ks 


(C.5) 


(C.6) 


—+ 3gl + —gj - {h^ + 6hl + hl)Sk3 - Xl^ 


X2 

'^kl2 


X2 

'^kl3 


-QiKlk + Kh + K]k + \lk) - n\li + \l + hi)} 


(C.7) 


+ 35'2 + ~ + ^6 + hl)^k3 - ‘^hldis 


X2 


X2 

^fc«2 


X2 

'^kl3 


-Q{\lk + \ll) - ^\]k - \lk - + ^2l + Asfc^z)} 


(C.8) 


^^'^‘^-^^kik — ^kik{-^9i + 3fi'2 + -^fl'3 ~ (2^6 + hl)5k3 - {hi + hl)5i3 - Xh^ - Xli2 - Xh^ 


«x'2 


'^'^kkl 


«X'2 

^'^klk 


9X'2 

‘^'^Ikk 


2(AJ, + A2,UA3^,)} 


(C.9) 


167r^^Aj^fc-A;^fc{-5f3 + 35'2 + ~ {hi + Qhl)Sk3 - hlSis 


X2 


X2 

^fc«2 


X2 

'^kl3 


fiX'2 

'^'^kkk 


X'2 


2a;^,, 


^\kk 


2(Aim + A2fcj + Ag^;)} - i/ \kkKi3h‘lhl5k3 (C.IO) 


A^^; — + 85 f| 2/i^((5fc3 + 2A^^^fc ‘^^kki ‘^^kik ‘^^ilk 


-6(a;^ + a;^ + a;'2j} 


(C.ll) 


167r^^^A2fc; — A 2 fci{ 3(71 + 85^1 2/i^((5fc3 + (5i3) ^X^n. ‘^^Xih. 


-Q{Kll + >^2ll + >^3ll)} 


(C.12) 


167r2^A3^i = A3|;{^^^ + 8c/| 


2/^? 


‘2hl{5k3 + 5«3) — ‘^^kkk ~ 2A 


2 

kkl 


9 x'2 

‘^'^klk 


9 x'2 
‘^'^Ikk 


-6(a;i, + a;^ + a;2,)} 


(C.13) 
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Figure 6: The running couplings of the RPV 13 system evaluated at the susy scale 

Msusy = ITeV (t ~ 61.25j, as a function of tanP, subject to the constraints mtimt) = 
165 GeV, mbimt) = 2.76 GeV, mrimt) = 1.78 GeV and unification of gauge couplings at 
Mgut = 2 X 10^® GeV. Figures a) and b) ( resp. c) and d) ) are two different configura¬ 
tions where the perturbativity bound is reached by one of the couplings at the GUT 
scale. 
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